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SUMMARY 

A direct solution procedure for computing the flutter Mach number and the 
flutter frequency is applied to the aeroelastic analysis of propfans using a 
finite element structural model and an unsteady aerodynamic model based on a 
three-dimensional subsonic compressible lifting surface theory. An approxi- 
mation to the Jacobian matrix that Improves the efficiency of the iterative 
process is presented. The Jacobian matrix is indirectly approximated from 
approximate derivatives of the flutter matrix. Examples are used to illus- 
trate the convergence properties. The direct solution procedure facilitates 
the automated flutter analysis in addition to contributing to the efficient 
use of computer time as well as the analyst's time. 


INTRODUCTION 

Flutter of propfans and other types of turbomachinery blading is an impor 
tant phenomenon that has generated considerable interest. Flutter prevention 
has been a significant factor in the design of propfan blades. Flutter preven 
tion is also significant for turbomachinery, particularly for unshrouded 
blades. As a result, a comprehensive flutter analysis procedure for rotating 
blade systems incorporating general mistuning and three dimensional subsonic 
aerodynamics is being developed at the NASA Lewis Research Center (refs. 1 
to 3). 

Research on the aeroelastic behavior of propfans has progressed from both 
the structural and the aerodynamic viewpoints. Because of the complex geome- 
try of the propfan, structural modeling has shifted from an initial beam model 
(ref. 4) to advanced finite element models (refs. 2, 3, 5, and 6). Aerody- 
namic models evolved from two-dimensional lifting line theory (refs. 2 and 7) 
to three-dimensional lifting surface theory (ref. 8). Development of aerody- 
namic models with supersonic, transonic, and separated flow effects is also a 
continuing effort at present. Aeroelastic analysis of propfans (refs. 2 to 5 
and 9) has kept pace with developments in structural and aerodynamic modeling. 
Comparisons with experimental data (refs. 2, 3, and 10) have resulted in rea- 
sonable confidence in analytical predictions. 
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With the recent advances in computer technology, automated design of prop- 
fan and turbomachinery blades using optimization techniques has become practi- 
cal. Efforts in this direction have already begun (refs. 11 to 15). Design 
optimization employing flutter constraints requires repeated solution of the 
aeroelastic equations of motion to obtain the flutter parameters as the design 
is updated. For the optimization to be performed in a realistic period of 
time, an automated flutter analysis capability is essential. It is also desir- 
able for the analysis to be computationally efficient in order to keep the cen- 
tral processing unit (CPU) time and the turnaround time within reasonable 
limits. 

The conventional flutter analysis procedure, based on the k-method 
(ref. 16), casts the equations of motion as a complex linear eigenvalue prob- 
lem where the real and imaginary parts of the eigenvalue are the effective 
damping and frequency respectively of the aeroelastic system. The solution 
consists of an inner-outer iteration loop to obtain the flutter Mach number and 
frequency. The inner iteration loop matches the assumed frequency to the cal- 
culated frequency while the outer loop iterates on the Mach number to obtain a 
zero effective damping. 

A stumbling block in the development of a reliable automated procedure has 
been the need to track several eigenvalues through possible crossovers. Previ- 
ous efforts in the development of automated flutter analysis procedures for 
fixed wings attempted to overcome this difficulty in two ways. One was to use 
eigenvalue derivative information to predict the crossovers and then keep track 
of them during the iteration (refs. 17 and 18). Apart from being cumbersome, 
this procedure required the derivatives of the aerodynamic force coefficients 
with respect to Mach number and assumed frequency. The other development was 
to use an accurate approximation to the eigenvalue which would retain the con- 
tinuity of the eigenvalues over the desired range of Mach number and assumed 
frequency (ref. 19). This latter method is susceptible to the accuracy and 
the reliability of the approximation over a wide range and requires careful 
monitoring. In addition, the inner-outer iteration loop in the conventional 
procedure is computationally expensive. The inner loop involves the computa- 
tion of the matrix of aerodynamic coefficients and the solution of a complex 
eigenvalue problem. This can degrade efficiency, particularly for problems in- 
volving structural and/or aerodynamic mistuning. 

The present work explores an efficient and automated procedure which 
replaces the costly inner-outer double iteration loop with a single iteration 
loop. The need to track the eigenvalues is also avoided. This is accom- 
plished by viewing the flutter problem as a complex implicit double eigenvalue 
problem and directly solving for the flutter Mach number and the flutter fre- 
quency by a quasi-Newton determinant iteration. The flutter eigenvalues are 
never calculated and the avoidance of eigenvalue tracking makes true automa- 
tion possible. An accurate approximation to the Jacobian matrix which makes 
the iterative process very efficient is a new feature of the present procedure. 
The proposed procedure is tested by applying it to the flutter analysis of a 
propfan model at a number of rotational speeds. 

While the idea of direct solution of the flutter equation is not new 
(refs. 16, 20, and 21), it has not been very popular. One reason may be that 
the alternative solution, involving flutter eigenvalues, has more physical 
appeal. However, the imperatives of automating the flutter analysis and reduc- 
ing the computational expense of an inner-outer iteration loop make the direct 
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solution more attractive. Some previous applications of direct solution of 
flutter equations by determinant iteration considered highly simplified aerody 
namic models. No results are available on the convergence and efficiency of 
the direct solution procedure considering more complex models such as kernel 
function aerodynamics, the computational expense of which makes such issues 
paramount. One of the objectives of this paper is to provide such results. 

The present procedure also Improves upon the previously published direct 
solution techniques of Hassig (ref. 16) and Cardani and Mantegazza (ref. 21). 
Hassig (ref. 16) uses a secant-based determinant iteration to obtain the 
matched flutter frequency and approximate subcritical rate of decay whereas 
the present procedure obtains both the matched flutter Mach number and the 
matched flutter frequency by a more efficient implementation of determinant 
iteration. Cardani and Mantegazza (ref. 21) solve the full set of nonlinear 
equations of motion without using determinant iteration to obtain the matched 
values of flutter Mach number, flutter frequency and the flutter mode simulta- 
neously. In the procedure developed here, only two nonlinear equations are 
solved to obtain the flutter point and the flutter mode is obtained very effi- 
ciently with the Information available from the final iteration. The Jacobian 
approximation and the direct solution procedure presented here are applicable 
to any aeroelastic stability analysis, including fixed wing, in a straightfor- 
ward manner. The Jacobian approximation is also generally applicable to deter 
minant Iteration in any application. 


NOMENCLATURE 

[A] generalized aerodynamic force matrix 

[B] system flutter matrix 

D system flutter determinant 

Dr real part of the flutter determinant 

Di imaginary part of the flutter determinant 

[J k ] Jacobian matrix 

k counter used in determinant iteration 

[Kg] generalized stiffness matrix 

M Mach number 

Mp flutter Mach number 

[Mg] generalized mass matrix 

AM k M P < k _i) - Mp < k) 

(qo> flutter mode of vibration 

(q> generalized coordinate vector 
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t time 


&0.75R 

blade setting angle at 

three-quarter radius 

Awk 

w F(k-l ) - “F(k> 


o r 

r-th inter-group phase 

angle 

Q 

rotational speed of the propfan 

0) 

frequency of vibration 


«F 

flutter frequency 



FORMULATION OF THE FLUTTER ANALYSIS PROBLEM 

The computational procedure to be presented is applied to the analytical 
formulation described in detail in references 2 and 3. This formulation is 
applicable to the flutter analysis of a single rotor propfan containing an 
arbitrary number of blades rotating at a fixed speed in an axial flow. The 
structure is modeled by finite elements. The aerodynamic model (ref. 8) is 
based on a three-dimensional subsonic compressible lifting surface theory. 

For simplicity, the effect of steady deformations due to aerodynamic loads 
on the flutter boundary is neglected. The error introduced by ignoring the 
steady aerodynamic deformations is studied in reference 2 and is shown to be 
small enough to produce an approximate flutter point. In particular, the 
effect of the steady aerodynamic deformations on the flutter Mach number is 
not very significant. Thus, the approximate flutter analysis neglecting steady 
aerodynamic deformations is suitable for use in design optimization procedures 
requiring repeated efficient execution of the flutter analysis. The optimal 
design can be easily checked for the flutter condition by using the refined 
flutter analysis with steady deformations and the conventional procedure. 

The propfan is assumed to have Identical groups of blades symmetrically 
distributed about a rigid disk. The linearized aeroelastic equations of 
motion are then uncoupled for different inter-group phase angle modes. The 
flutter Mach number for the propfan is then the lowest Mach number at which 
one of the intergroup phase angle modes becomes unstable. For a given inter- 
group phase angle a r , the equations of motion, after modal transformation, 
can be written in the form (ref. 3): 

[M g ]{q} + C Kg ] {q} = CA(M,to)]{q} (1) 

where [Mg] is the generalized mass matrix, [Kg] is the generalized stiff- 
ness matrix, {q} the generalized coordinate vector and [A(M,u>] the general- 
ized aerodynamic matrix for a group of blades. The aerodynamic matrix 
[A(M,u> ] is usually valid only for simple harmonic motion of the airfoil. For 
full details of the construction of these quantities, the reader is referred 
to references 2 and 3. 
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For simple harmonic motion. 


{q> - {q 0 >e'“ t 


( 2 ) 


Equation (1) now leads to 

<-<o 2 CM g ] + [K g ] - [A(M,w)]){q 0 } = {0} (3) 

The conventional procedure for obtaining the flutter Mach number and fre- 
quency, outlined In references 2 and 3 and Illustrated In figure 1 (reproduced 
from ref. 2 Ignoring the computation of steady deformations due to aerodynamic 
loads). Is as follows. The aerodynamic matrix Is evaluated at an assumed Mach 
number and an assumed frequency and then equation (3) Is solved for all the 
eigenvalues, co 2 . This procedure has been Implemented for propfans In a pro- 
gram called ASTR0P3. In general, these eigenvalues are complex. The real and 
Imaginary parts of ico represent the effective damping and frequency, respec- 
tively. The assumed frequency Is varied until It Is equal to the frequency 
corresponding to the eigenvalue with the least effective damping. This fre- 
quency matching forms the Inner Iteration. When this Iteration reaches conver- 
gence, the Mach number Is varied until the effective damping of the eigenvalue 
corresponding to the matched frequency Is equal to zero. This forms the outer 
Iteration. The flutter Mach number and the flutter frequency are obtained at 
the convergence of the outer Iteration. 

The conventional procedure cannot be reliably automated because It 
requires the preservation of the Identity or ordering of the eigenvalues over 
a wide range of assumed frequencies and Mach numbers. Most eigensolution rou- 
tines do not compute the eigenvalues in any particular order and the sorting 
of eigenvalues by frequency or magnitude does not usually preserve the continu- 
ity. Loss of continuity necessitates user intervention and complicates the 
automated analysis. 

In order to alleviate these problems, a direct solution of equation (3) 
viewed as an implicit double eigenvalue problem Is proposed and described 
below. 


DIRECT SOLUTION OF THE FLUTTER EQUATION 

When the dependence of the aerodynamic matrix on the assumed Mach number 
and frequency Is considered explicitly, equation (3) can be viewed as an 
implicit double eigenvalue problem. If Mp is the flutter Mach number and 
cop the flutter frequency, then equation (3) leads to 

[B]{q 0 } = {0} (4) 

where CB] = -co F 2 [M g ] + CK g ] - [A<Mp,<op>] 

In general, the aerodynamic matrix [A(Mp,<op>] Is a transcendental func- 
tion of the assumed frequency and Mach number. Mp and cop are the two eigen- 
values of equation (4). Only real values of Mp and cop are of Interest. 

The two eigenvalues Mp and cop are coupled since the aerodynamic matrix is 
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complex. Cardani and Mantegazza (ref. 21) proposed a direct solution of equa- 
tion (4) and a normalizing condition, treating Mp, up and the components of 
{qO} as unknowns. For a nontrivial flutter mode, we have the characteristic 
equation 

det (-up 2 CMg] + [K g ] - [A(Mp.up)]) = 0 (5) 

Let 

D = det ( —up 2 [M g ] + [K g ] - [A(Mp,u F )]> 

= Dp(Mp,up) + iDi<Mp,up) (6) 

where Dr and Dj are the real and the imaginary parts of the characteristic 
determinant D respectively. Then, at flutter condition. 


DR(Mp,up) = 0 
Dj(Mp.up) = 0 


(7) 


Equations (7) consist of two equations in two unknowns, Mp and up. 

These equations can be solved by any of the methods for solving a system of 
nonlinear equations. When equations (7) are solved for Mp and up, equa- 
tion (4) is satisfied and no further iterations are required for the purpose 
of matching assumed and computed quantities. This procedure is illustrated in 
figure 2. Once Mp and up are found, inverse iteration can be used to find 
the flutter mode. 


In contrast to the conventional procedure, the direct solution outlined 
here eliminates the need to track eigenvalues to determine the flutter point. 

In addition to this important benefit, the double iteration on the complex 
eigensolution is replaced by a single solution of a system of two real nonlin- 
ear equations. The Mach number and frequency are varied simultaneously in 
this procedure rather than one at a time as in the conventional procedure. 

Thus the flutter Mach number and the flutter frequency are determined 
simultaneously. 

The formulation of a transcendental double eigenvalue problem in prefer- 
ence to a linear single eigenvalue problem may seem to defeat the objective of 
increased efficiency, even if it is more suitable for automation. However, 
the price to be paid is not as great as it may seem because the transcendental 
eigenvalue problem needs to be solved for only one set of eigenvalues in most 
cases, whereas the linear eigenvalue problem has to be repeatedly solved for 
all the eigenvalues which equal in number to the number of assumed mode shapes. 


NUMERICAL SCHEME FOR DIRECT SOLUTION 

Equations (7) may be solved for Mp and cop by Newton's method. The 
iterative scheme for Newton's method is 


k | 

) l 

= 

( Mp) 

- W 1 

1 DR j 

f W F ' 

) (k+l ) 1 

U F \ 


[ Dl 1 

1 


1 1 (k) 
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where k is the iteration number and [J^] is the Jacobian matrix and is 
given by 


[J k ] 


kR,M f 

^kl ,Mr 


'kR,» F 

^I.Wr 


( 9 ) 


The Jacobian [J^] is expensive to compute because the evaluation of the 
aerodynamic matrix [A^] is computationally intensive. Several quasi-Newton 
algorithms which approximate the Jacobian in various ways are available. In 
the following, a quasi-Newton algorithm is proposed that is more efficient than 
the currently available algorithms for the determination of the flutter Mach 
number and the flutter frequency. The numerical scheme is based on the hypoth- 
esis that approximating the derivatives of the flutter matrix [B^l provides a 
more accurate approximation to the Jacobian matrix [J|<] than directly approx- 
imating the derivatives of the characteristic determinant. The numerical 
scheme based on this hypothesis approaches the Newton's method in its superior 
convergence characteristics with the same cost per iteration as the secant 
method . 


The iterative scheme of equation (8) is first transformed in order to 
explicitly determine the relation between the derivatives of the characteris- 
tic determinant and those of the flutter matrix. This is accomplished by 
employing the following remarkable result, known as the Trace theorem, proved 
by Lancaster (ref. 22 ): If the elements of C B( a) ] are differentiable func- 

tions of a, then for any a for which the characteristic determinant 
D(a) * 0 , 


0,5 = D • trace of ([B] _1 tB],*) 

When a Is successively set to Mp and up, equation ( 10 ) leads to 
£>k,Mp - D kR,M F + iD kI,M F * D k * trace of ([B k ]- ] CB k ] >Mp ) 
Dk,o, F - D kR,«p + 1D kI,wp = • trace of ([B k H CB k ],„ F ) 

The Iterative scheme of equation (8) can then be written as 


M, 


»r 


(k+1 ) 



1 

( 9r V ! '“F " 9R ’“F 9l ’ M 0 (k) 


9 I,«r 

- g I,M F 



where 


(g.Mp) - trace of ([B^] -1 , H p> 

f9»«p') - trace of ([B^H [B^l , w p) 
v J k 





(10) 


( 11 ) 


( 12 ) 


( 13 ) 
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The Iterative formula of equation (12) does not require the determinant of 
the coefficient matrix [B^] of the flutter equation (4). However, is 

required to compute g.Mp and g, U p. The computational cost for the inversion 
of [Bfc] is usually negligible in comparison to that of forming the matrix 
C8 k l. 


The derivatives CB k D ,Mp and CB|<], UF are approximated by following a 
reasoning similar to that employed in deriving Broyden's method (ref. 23). 

Let AMk = Mp(|<_D - Mp(k) and Aw|< = wF(k-l) - «F(k)* The derivatives are 
approximated in the direction of the last move to satisfy 

C B k _ i ] = [ B k ] + CB k ] ,Mp * AM k + C B kl»up * ^k (14) 

and are assumed to be unchanged in the direction orthogonal to the last move 
and this Implies 

[Bk],Mp<MF - Mp(|<.)) + [Bk], UF ( “F ~ “F(k)> = C B k- 1 ^ ,M F <M F - M F(k) > 

+ [Bk-l ^uf ( “F - “F(k)> <15) 


for any Mp and top satisfying: 

(Mp - Mp ( k ) ) • AM^ + (top - w F(k)) * A“k = 0 
Equations (14) and (15) are solved to obtain 

UB k _,] - [B k ]>AM k ♦ <CB k _,3. MF • • 4\)A- k 


[B,], 


k ’M, 


[B kV 


AMj^ + Acoj; 


(16) 


(CB k -l3 - [B k ]) H * <[B k-l 3 *M F ‘ H - [B k-l ] *<o F * AM k )AM k 


AMj^ + Auj^ 


(17) 


It was expected that equations (13), in conjunction with equations (16) 
and (17), approximate the Jacobian matrix accurately enough to obviate the 
need to update the Jacobian in the entire range of convergence. Results in 
the next section show that this is Indeed the case for the examples studied, 
even when large moves were made during the iteration. 

The proposed numerical scheme for direct solution of the flutter equation 
is: 


(1) Obtain Initial estimates Mp(o) and uf( 0) f° r ttie flutter Mach 
number and flutter frequency. 

(2) Obtain initial estimates for the derivatives [Bo],m f and ^O^'uf- 
These are best obtained by computing finite difference approximations. 

(3) Compute new estimates Mf(k+ 1 ) and w F(k+D* k = 0, 1. 2,... using 
equation (12) . 
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(4) Test for convergence and stop if converged. 

(5) Increment k * k + 1 . 

(6) Compute updated derivatives C6^3 .Mp and [B^],^ using 
equations <16> and (17). 

(7) Repeat from step 3. 

[Bfc]- 1 from the final Iteration of the above procedure can be used In an 
Inverse Iteration scheme to determine the flutter mode at little additional 
cost. Thus, If {qop) Is the required flutter mode and {v} Is any vector 
that Is not totally deficient In {qop). then the following Inverse Iteration 
scheme Is used to obtain the flutter mode: 

{qOF>l * CB] _1 {v} 

(QOFJm+l ■ CB3 -1 (dOF^m* m - 1. 2,... (18) 

Normally, one step of the Inverse Iteration scheme Is sufficient to obtain 
flutter mode within a reasonable accuracy. 


RESULTS AND DISCUSSION 


To demonstrate the current procedure, flutter boundary calculations were 
performed for two configurations of propfan rotors described In references 2 
and 3. The first configuration consists of eight Identical blades, designated 
SR3C-X2 , made of graphl te-ply/epoxy-matrl x material. The second configuration 
Is an alternately mlstuned rotor with eight blades. The rotor was constructed 
by alternately mounting two sets of composite blades. The first set of four 

blades, SR3C-X2, are the same as those on the rotor of the first configuration. 

The second set of blades, SR3C-3, are Identical to the first set except for the 
ply angles of the composite material. For details on the blade models, the 
reader is referred to references 2 and 3. In this report, the first configura- 
tion is called the tuned rotor and the second configuration the mlstuned rotor. 

To validate the current procedure and the computer program, flutter bound- 
aries were calculated at various rotational speeds for both the tuned and mis- 
tuned rotors. Convergence was considered to have been obtained when the 
solution did not change to four significant digits In successive iterations. 
During the Iteration process, the value of the Mach number sometimes turned 
out to be supersonic at the tip. No aerodynamic analysis can be performed at 
such Mach numbers using the current aerodynamic model. When such Intermediate 
solutions (including those Involving negative Mach numbers) arose, the move 
length was cut in steps of 80 percent until an acceptable Mach number was 
obtained. This technique was found to be effective In driving the solution to 
the desired range without an excessive penalty in the number of Iterations. 

The results obtained with the present procedure starting with various Initial 
guesses are compared with the "exact" flutter boundaries obtained using the 
conventional procedure, at two rotational speeds each for tuned and mistuned 
rotors. The flutter boundaries are in excellent agreement for these cases, 
thus validating the accuracy of the procedure and the convergence criteria. 
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The typical progress of iteration, for the initial guesses for flutter 
Mach number and flutter frequency of 0.5 and 310 Hz respectively, with the 
present direct solution procedure and the conventional procedure, is shown in 
table I. Note that the conventional procedure relies on user interaction and 
judgment. The progress of iteration shown for the conventional procedure is 
typical. The current procedure, in addition to being suitable for automation, 
is also more efficient as evidenced by the considerably smaller number of 
analysis steps. Thus both the CPU time and analyst man-hours are considerably 
reduced using the present procedure. 

To judge the accuracy of the approximate Jacobian, the nonlinear equa- 
tions (7) were solved by Newton's method, the present numerical procedure and 
by alternate quasi-Newton methods such as the multi-point secant method 
(ref. 24) (Implemented in IMSL routine ZSCNT), modified Powell's algorithm 
(Implemented in IMSL routine ZSPOW) and by Broyden's method (ref. 23). For the 
cases tested, the present procedure outperformed all three alternate methods 
in terms of efficiency. Figures 3, 4, and 5 illustrate the convergence behav- 
ior of the current procedure in comparison to Newton's method and the multi- 
point secant method. For the cases shown, the "exact" flutter Mach number was 
0.641 and the "exact" flutter frequency was 294 Hz. Even though the character- 
istic determinant is never calculated in the current procedure, the varia- 
tion of the absolute value of with each Iteration is shown in figure 3 for 
comparison with the Newton's method and the secant method. The Initial values 
for Mp and top were 0.65 and 330 Hz respectively. The determinant value has 
been scaled so that 1.0 < |Dq| £ 10.0. It is seen that the iteration history 
using the current numerical procedure closely resembles that of the Newton's 
method indicating the accuracy of the approximation proposed here for the 
Jacobian matrix. In contrast, the secant approximation for the Jacobian 
matrix almost doubles the number of iterations required for convergence in 
this case. 

The range of convergence is an important factor in any Iterative proce- 
dure since it has an Important effect on how closely the initial solution must 
approximate the final solution. Figure 4 shows the number of flutter matrix 
evaluations required for convergence for various initial guesses for Mp and 
a fixed initial guess for «p = 310 Hz. The range of Mach number convergence 
is seen to be from Mp = 0.2 to Mp = 0.8. Figure 5 similarly shows the 
number of flutter matrix evaluations required for convergence for a fixed ini- 
tial guess for Mp of 0.65 and various Initial guesses for top. The frequency 
range of convergence with the current procedure, top = 230 Hz to top = 350 Hz, 
is seen to be slightly larger than that with the secant method. From these 
figures, it can be stated that the present procedure has a large radius of 
convergence. 

The results obtained show that a fair initial guess would converge to the 
"exact" flutter point at the expense of about 5 to 10 flutter matrix evalua- 
tions. Table II shows the CPU times on the CRAY-XMP that were required to 
obtain the flutter boundary for the case of a good initial guess and a poor 
initial guess. The CPU times for one flutter eigenvalue analysis at a given 
set of Mach number and assumed frequency are also shown for comparison. With 
a good initial guess, the flutter Mach number and the flutter frequency can be 
obtained for two or three times the cost of a single eigenanalysi s . The 
present procedure is much less expensive in terms of CPU time as well as ana- 
lyst's man-hours than the conventional procedure, even though precise compari- 
sons have not been made. However, the primary motivation for the current 
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procedure has been the feasibility of automation rather than the computational 
expense. 

The current procedure may not find the lowest flutter Mach number, if more 
than one structural modes were to flutter in the Mach number and frequency 
range of interest for the selected intergroup phase angle mode. Under these 
circumstances, one will be forced to search the entire range of interest for 
the roots of the equations (7) starting with different initial guesses and this 
is not amenable to an efficient automated procedure. However, it is expected 
that such is rarely the case for tuned or alternately mistuned propfans. This 
is not a major limitation for two other reasons: (1) the frequency interval in 

which flutter speeds occur is usually determined early in the design phase and 
(2) the search domain can be considerably reduced after a few orientational 
runs. 


CONCLUSIONS 

A direct solution of the equations of motion is demonstrated to be a reli- 
able automated flutter analysis procedure if steady aerodynamic deformations 
are ignored. The direct solution method replaces the inner-outer iteration 
loop of the conventional procedure by a single iteration loop. A numerical 
procedure, based on an accurate and efficient approximation to the Jacobian 
matrix, is presented. The procedure is straight-forward in concept and results 
for test cases show good convergence properties. Since the procedure is itera- 
tive, it is particularly suitable for design optimization. This is because, 
as the optimal design is evolved, the flutter solution is expected to change 
incrementally from design to design, so that the previous solution provides 
good estimates for the current solution. 
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TABLE I. - PROGRESS OF ITERATION 
[£2 = 5280 rpm, 3q.75R - 61.6°, o r = 225°.] 


Count 

Conventional procedure 

Present procedure 

M 

G>, 

Hz 

M 

to, 

Hz 

1 

0.500 

310.0 

0.500 

310.0 


.500 

267.5 

.499 

310.0 


.500 

268.9 

.500 

313.0 


.500 

268.9 

.701 

289.7 

5 

.700 

268.9 

.590 

287.7 

6 

.700 

299.8 

.641 

293.6 

7 

.700 

298.9 

.642 

294.1 

8 

.700 

298.9 

.641 

294.1 

mm ■ 

.616 

286.3 

a . 641 

294.1 

mm 

.616 

290.1 



in 

.616 

290.5 



ns 

.616 

290.4 




.640 

292.9 



Mm 

.640 

293.9 



B9 

.640 

294.0 



IH 

.641 

293.9 



ul 

.641 

294.1 



mm 

a . 641 

294.1 




Converged. 


TABLE II. - CPU SECONDS FOR AUTOMATED FLUTTER ANALYSIS 



Eigensol uti on 
at a single 
set of 

Direct solution to find 
Mp and cop 


M and u 

Good 
ini tial 
guess 3 

Poor 
i ni ti al 
guess 3 

Tuned rotor b - 8 blades; 

6 modes/blade; 5280 rpm; 
o r = 225° 

4.332 

10.356 
(M = 0.70 
co = 310 Hz) 

22.146 
(M = 0.45 
co = 340 Hz) 

Mistuned rotor c - 8 blades; 
4 groups; 2 modes/blade; 
5190 rpm; o r = 90° 

10.020 

22.084 
(M = 0.65 
co □ 310 Hz) 

31.970 
<M = 0.5 
co = 340 Hz) 


^Initial guesses are given in parentheses. 
b "Exact" Mr = 0.641 and "exact" cop = 294 Hz. 
c "Exact" M F = 0.718 and "exact" « F = 285 Hz. 
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CHANGE M, 




FIGURE 2. - AUTOMATED PROPFAN FLUTTER ANALYSIS. 




ABSOLUTE VALUE OF THE 
FLUTTER DETERMINANT, |D| 



ITERATION NUMBER 


FIGURE 3. - VARIATION OF THE ABSOLUTE VALUE 
OF THE DETERMINANT WITH ITERATION NUMBER 
(H - 5280 RPM, f} 0 .75R* 61.6°, 0 r =225°). 



FIGURE 4. - CONVERGENCE OF DIRECT SOLUTION 
FOR VARIOUS INITIAL GUESSES AT FLUTTER 
MACH NUMBER. 



ASSUMED FLUTTER FREQUENCY, Up(Hz) 

FIGURE 5. - CONVERGENCE OF DIRECT SOLUTION 
FOR VARIOUS INITIAL GUESSES AT FLUTTER 
FREQUENCY. 
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